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Abstract. We note that two formulations of dual gonihedric Ising models in 3d, one 
based on using Wegner's general framework for duality to construct a dual Hamiltonian 
for codimension one surfaces, the other on constructing a dual Hamiltonian for two- 
dimensional surfaces, are related by a variant of the standard decoration/iteration 
transformation. 

The dual Hamiltonian for two-dimensional surfaces contains a mixture of link and 
vertex spins and as a consequence possesses a gauge invariance which is inherited by 
the codimension one surface Hamiltonian. This gauge invariance ensures the latter is 
equivalent to a third formulation, an anisotropic Ashkin-Teller model. We describe the 
equivalences in detail and discuss some Monte-Carlo simulations which support these 
observations. 
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1. Introduction 

The dual of the standard Ising Hamiltonian with nearest neighbour (ij) couplings on a 
3d cubic lattice 

H ising = — cr i<3 j (1) 

(ij) 

is the Z 2 Ising gauge theory 

Hcauge = _ UijUjkUklUu (2) 

[ijki] 

where the sum is over plaquettes [ijkl] and the spins live on the edges of the lattice. 
The coupling (3 in the partition function Z{0) = 2\^{o-} ex P( — fiHi sing ) and its dual (3* in 
Z (P*) = E {U } e M-P*H Gauge ) are related by 0* = -(1/2) log tanh/3. 

In this paper we will investigate the relation between three (apparently) different 
formulations of the dual to the gonihedric Ising model [l] 

H k= q = — OiOjO k oi (3) 

[ijki] 

which, like the Ising gauge theory, has a plaquette interaction but in which the spins 
now reside at the vertices of the 3d lattice. The subscript k = appears because 
this plaquette Hamiltonian is a particular case of a one-parameter family of gonihedric 
Hamiltonians [|] 

Hgonihedric = ~4k ^ (JjO j + K ^ (J iO j — (1 — k) OiOjO k Oi . (4) 

(ij) ((ij)) [ijkl] 

defined by Savvidy and Wegner [2], where the ((ij)) are next-to- nearest neighbour 
sums. The spin cluster boundaries of this Hamiltonian were intended to mimic a 
gas of worldsheets arising from a gonihedric string action. When discretized using 
triangulations, this action may be written as 

5 = ^1^-^10(0^), (5) 

where 6(a>ij) — 1 7T Otj^j | - Ot<ij IS the dihedral angle between the neighbouring triangles 
with a common edge (ij) and \X{ — Xj\ are the lengths of the triangle edges. 

The word gonihedric was originally coined to reflect the properties of this action 
which weights edge lengths between non-coplanar triangles rather than their areas. It 
combines the Greek words gonia for angle, referring to the dihedral angle, and hedra for 
base or face, referring to the adjacent triangles. H gonihedric is an appropriate cubic lattice 
discretization of such an action because it too assigns zero weight to the areas of spin 
cluster boundaries, rather weighting edges and intersections [3j. This gives H gonihedric 
very different properties to Hi S i ng where (only) the areas of spin cluster boundaries are 
weighted. 



| We have dropped a factor of 1/2 in the coupling definition compared with [5|[6] in order to keep the 
standard definition of the duality relations here. 
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The plaquette action H K=0 has been shown to possess a degenerate low-temperature 
phase and a first order phase transition as well as interesting, possibly glassy, dynamical 
properties [4]. It displays a peculiar "semi-global" symmetry in which planes of spins 
may be flipped at zero energy cost, accounting for the degeneracy of the low temperature 
phase. For non-zero K this symmetry appears to be broken at finite temperature and 
the transition becomes second order. In |5] we observed that one formulation of the 
dual to H K=0 , which took the form of an anisotropic Ashkin- Teller model, displayed 
similar symmetry properties since it was possible to flip planes of spins in this also. 
In [6] we related this to a superficially different dual formulation that employed three 
flavours of spins by using a gauge-fixing procedure. There were indications of potentially 
interesting dynamical behaviour for both Hamiltonians in [5] and [6]. 

In the following we discuss the derivation of these two dual Hamiltonians as well 
as a third possibility. We then describe the relation between the various Hamiltonians 
via a decoration transformation and gauge-fixing, before outlining some Monte-Carlo 
simulations in support of these observations. 

2. Duals Galore 

The dual to H K=0 was initially constructed "by hand" by Savvidy et.al. [7] by considering 
the high temperature expansion of the plaquette Hamiltonian 

Z((3)= J2 e M-(3H K=0 ) 

M 

= yi n c ° sh c 9 ) t 1 + tanh c 9 ) (wjwi)] ( e ) 

{a} [ijkl] 

which can be written as 

Z{0) = [2 cosh (/3)] 3L3 ^nh (/3)] n(S) (7) 

{5} 

on an L 3 cubic lattice, where the sum runs over closed surfaces with an even number of 
plaquettes at any vertex. In the summation n(S) is the number of plaquettes in a given 
surface. 

Surprisingly, the low temperature expansion (i.e. high temperature in the dual 
variable /3* = —(1/2) log tanh 0) of the following anisotropic Hamiltonian 

HdualO = —J"] ^(Tj - ^ T i T k ~ X] W ( 8 ) 

(ij) (ik) (jk) 

produced the requisite diagrams. In Hd ua io the sums are one- dimensional and run along 
the orthogonal axes, with ij, ik and jk representing the z, y and x axes respectively 
using our conventions. The spins are non-standard and live in the fourth order Abelian 
group, since the geometric constraint on the plaquettes means that 

ea = a , er = r , er\ — rj 

a 2 =r 2 =r] 2 = e (9) 

<7T = Tj , TTj = (X , T]0~ — T 
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with e being the identity element. They can be thought of as representing differently 
oriented matchbox surfaces such as that shown in Fig. (f), which are combined by 
facewise multiplication. The shaded faces carry a negative sign and the associated spin 
variable lives at the centre of the matchbox. Any spin cluster boundary in the model 
can be constructed from such matchboxes while still satisfying the local constraint on 
the number of incident plaquettes. 




Figure 1. An elementary matchbox surface represented by one of the spins in equ. K\ 



The spins may also be taken to be Ising (±1) variables if we set r]i = <Ji Ti, which 
is more convenient for simulations. This modifies Hd ua io to an anisotropically coupled 
Ashkin- Teller Hamiltonian [8] 

Hduall = - 2j a i°j ~ ^ TiTk ~ S ^ a i akT i Tk ■ (10) 

(ij) (ik) (jk) 

We investigated this formulation of the dual model in [5] and found that it displayed a 
first order phase transition and similar semi-global symmetries to those of H k= q. The 
symmetries were a direct consequence of the anisotropic couplings, which allowed a 
greater freedom in transforming the spin variables than in the isotropically coupled 



version of equ. (10), which is just the Ashkin- Teller model at its four-state Potts point. 

It is also possible to construct duals to H K=0 and its higher dimensional equivalents 
[9] using the completely general framework for duality in lattice spin models that was 



originally formulated by Wegner in 10 . There are two possible ways to write the dual 
to H K= o in three dimensions with this machinery, using either the general formula for the 
dual of codimension one surfaces or the formula for the dual of two dimensional surfaces 
in d dimensions. If we temporarily use the notation of [9], the dual Hamiltonian for a 
codimension one surface in d dimensions is given by 

H dual ,codiml 

a<(3,r 7 

(11) 

where the A spins live on each of the (d — 3) dimensional (hyper) vertices situated at 
the vertices r of the hypercubic lattice and the indices a, (3, 7 run from 1 to d. The unit 
vectors e 7 point along the lattice axes. The dual Hamiltonian for a two-dimensional 
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gonihedric surface embedded in d dimensions is of the form 

HLi,2d = -X)X>07(Wr,f+e 7 )A0 7 (f+e 7 ) (12) 

f ^7 

where we now have T spins on each (hyper)edge in addition to the A spins at each 
vertex. 

If we specialize to two dimensional surfaces embedded in three dimensions, which 
is the case for the dual of H k= q, either formulation may be employed. Returning to 



our own notation [6], the codimension one Hamiltonian of equ. (11) in three dimensions 
may be written as 

H d ual2 — — 2_j ViVjUiHj - ^ T i T k^i^k ~ / ] <?j<?kTjT k , (13) 

(ij) (ik) (jk) 

where we again have one-dimensional sums as with Hduaio and Hd ua ii, but there are now 
three flavours of spins living at each vertex which display a local Ising gauge symmetry 
<jj, Tj, Hi — > ji/Ji, jiTi, jifa in addition to the planar flip symmetries of H dua ii- 

Still within the general approach of Wegner 10 , in three dimensions the 



Hamiltonian of equ. (12) for the two-dimensional surface variant also contains three 



flavours of vertex spins <7j, Tj, /ij, but in addition there are spin variables L^' 2 ' 3 living on 
the lattice edges which couple in an anisotropic manner to the vertex spins 

HdualS = ~Y1 ( a i U ij a 3 + ViU?jVj) ~ Yl ( T i U ik T k + Vi U ikVk) 

(ij) (ik) 

~ >I(M>a- • r ; r;U) • (14) 

(ifc) 

We thus have three superficially rather different Hamiltonian formulations for the, 
presumably unique, dual of the plaquette Hamiltonian H K=0 in three dimensions: 



Hduai3 in equ. ( 14 ) containing both vertex and edge spins 



Hduai2 m equ. (13) containing purely four spin interactions 



Hduaii in equ. (10) which is Ashkin- Teller in form. 



In the the next section we discuss the relation between Hd ua i3 and Hd ua i2, and thereafter 
that between H dua i2 and H dua ii- 



3. Decoration 

The equivalence between H dua i3 and H dua i2 is a consequence of a variation of the classical 
decoration transformation [9 11. In the standard transformation an edge with spins 
o"i,(T2 at each vertex is decorated with a link spin s as in Fig. (2). If the coupling 
between s and <j\ and 02 is (3, summing over the central spin s gives rise to a new 
effective coupling between the primary vertex spins a\, 02 



E 



cxp 



/3s(ai + a 2 ) = A exp(/3a 1 a 2 )- 



(15) 
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Figure 2. The standard decoration transformation 



Both the prefactor A and the coupling may be expressed in terms of ft by enumerating 



possible spin configurations in equ. (15). This gives 



A = 2cosh(2/3) 



1/2 



/3 = -logcosh(2/3). 



(16) 



/3 (<J 1 Ul 2 a 2 + HiUl 2 n 2 ) = A exp(pcri<j 2 fJ>ilJ>2)- 



(17) 



We can repeat this procedure with the U spins on each edge in H dua i 3 . In this case each 
direction has two flavours of vertex spin and performing the sum generates the four-spin 
couplings of Hduai2, for example 

exp 

The sum over U may be carried out globally over every edge which immediately 
demonstrates equivalence of the partition functions for H dua i 3 and H dua i 2 

Z = ^2 exp[-(3H dual3 \ 

{U,a} 



£ 

{U,a} 



cxp 



(ij) (ik) 



= B Y eX P P ^ WjtkPj + T iTkfUHk + Y °i°k T 3 T k 

M \_ (ik) (jk) 

= B ^2exp[-(3H dual2 \. 

M 

The overall factor B coming from a product of A's on the individual links is irrelevant for 
calculating physical quantities and the two couplings are again related by the standard 
decoration relation, (3 = (1/2) logcosh(2/3). 
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4. Gauge Fixing (and Flips) 

The equivalence between H dua i 2 and H dua n, on the other hand, is a consequence of the 
additional gauge symmetry [6j which is present in Hd ua i2 

a h n, Hi ->■ 7iCii, 7^^, 7i/ii . (19) 

We are at liberty to choose the Ising spin gauge transformation parameter ji to be equal 
to one of the spin values, say /ij, at each site so the gauge transformation then becomes 

o-j, n, ^ -> Wi, Hi T u 1 (20) 

which, using the fact that the sum over the remaining spin variables <jj, r% is invariant 
under the transformation, relates the partition functions for the two Hamiltonians as 

Z = exp {-(3H dual2 (a, r, //)] 

{(T,T,tl} 

= 2 L " ex P l-PH dual2 (a, t,/i= 1)] (21) 

W,r} 

= 2 L " ex P [-pH duall (a, r)] . 

The coupling /3 is not transformed in this case and we can, of course, choose to eliminate 
any one of the three spins, which simply amounts to relabelling the axes. From this 
perspective H dua n is simply a gauge-fixed version of H dua i 2 . This can be confirmed by 
Monte-Carlo simulations which measure the same energies (and energy distributions) 
and transition points for the observed first order phase transitions [6]. 

The equivalence between H dua i 3 and H dua i 2 described in the preceding section via 
the decoration transformation also sheds light on the somewhat unexpected presence of 
this gauge symmetry in H duai2 . All the terms in H dua i 3 are of the gauge-matter coupling 
form aiUijUj, so this action possesses a similar, standard gauge invariance to that seen 
in other gauge-matter systems such as the Z 2 gauge-Higgs model, namely 

rr!,3 . rrl,3 

r, -)• 7.T, , Tj -)• 7,-r,- , f/J 3 -> 7i^3' 3 7i (22) 

When the C/ spins are summed over to give H dua i 2 , the gauge symmetry of the a, r and 
spins remains as an echo of this symmetry. In both cases if we look at a single site 
transformation all three spins (Tj, Tj and Hi must be transformed. In H dua i 3 this is a 

12 3 

consequence of the way in which the three edge spins ' couple to the vertex spins. 

A characteristic feature of both H dua i 2 and H dua n is the flip symmetry of the low 
temperature phase, which allows planes of spins to be flipped at zero energy cost. This 
can be observed in the ground state by decomposing the full lattice Hamiltonian into 
cube terms and searching for minimum energy configurations on a single cube. The 
absence of any sign of non-zero magnetic order parameters in the low temperature phase 
in Monte-Carlo simulations then indicates that this ground-state symmetry persists 
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throughout the low temperature phase. More specifically, for Hd ua i2 it is possible to flip 
planes of pairs of spins at no energy cost. Similarly, planes of either one or two spins 
depending on the orientation may be flipped in Hduaix- Similar behaviour is also seen 
with the original H k= q plaquette Hamiltonian, where the symmetry has been confirmed 
to persist into the low-temperature phase in low-temperature expansions by Pietig and 
Wegner [12). 

In H dua i 3 flipping the three spins at a vertex may be compensated locally by flipping 
the six incident edge spins as shown in Fig. (3), which is the local gauge transformation 
that is still present in H dua i 2 . On the other hand, if just two spins are flipped at the 

+++ 

+++ 



+++ 




+++ 



+++ 



+++ 

Figure 3. Flipping three spins at a site may be compensated by flipping the six 
incident edge spins shown in bold, a purely local gauge transformation. 



central site, e.g. a and /x, the gauge transformation can no longer be applied to keep 
the disturbance local. However a global, planar spin flip can still leave the energy 
unchanged, as shown in Fig. (4). Choosing to flip different pairs of spins at the central 
vertex can be compensated by flipping the appropriate pairs of incident edge spins {U 2 
or U 3 ) and differently oriented planes of vertex spins. The planar flip symmetry is thus 
still a feature of the ground state of H dua i 3 and distinct from the gauge symmetry. 



5. Monte Carlo 



Monte-Carlo simulations reveal a first order phase transition in both H dua i 2 and H dua n 
at P ~ 1.39 [5}[6]. This can be seen in measurements of the energy, where there is a sharp 
drop at the transition point. A plot of the energy is shown for various lattice sizes in 
Fig. ([5]) for H d uai2, the values for Hd ua ii are essentially identical. The first order nature 
of the transition for H du ai2 and Hduaix is confirmed by observing a dual peak structure in 
the energy histogram P(E) near the transition point and a non-trivial value of Binder's 
energy cumulant 

(E A ) 

"»= 1-3^1 (23) 
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Figure 4. Flipping two spins (<r, // in this case) at a site may be compensated by 
flipping two incident edge spins (U 1 ) shown in bold along with the other coplanar 
vertex spins. This is no longer a purely local gauge transformation since the motif 
must be propagated across the lattice to maintain the same energy. 



-0.8 




-1.5 1 1 1 ' ' 1 1 1 

1.26 1.28 1.3 1.32 1.34 1.36 1.38 1.4 

Figure 5. The energy for Hduaii on lattices ranging from 10 3 to 20 3 from left to right. 
The lines joining the data points are drawn to guide the eye. Data from Hduaii is 
essentially identical. 

as a consequence of the shape of P(E). 

Based on these observations, and allowing for a factor of 1/2 in our definitions of 
Hduaii and Hduan in (5j[6] , we would expect to see a transition in Hduaiz at the the value of 
(3 found by inverting the decoration transformation, namely (1/2) ■ cosh _1 (exp(1.39)) = 
1.034 in the thermodynamic limit. To confirm this expectation, we carried out Monte- 
Carlo simulations using 10 3 , 12 3 , 16 3 and 18 3 lattices with periodic boundary conditions 
for all spins at various temperatures with a simple Metropolis update. After an 
appropriate number of thermalization sweeps, 10 7 measurement sweeps were carried 
out at each lattice size for each temperature. We have not attempted to construct 
a cluster algorithm since they do not offer effective speedup at first order transition 
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points such as that (presumably) under investigation here and because of the additional 
complication of the edge spins. 

Looking at measurements of the energy from our simulations of Hduaiz in Fig- ^ 
we can see that a similar sharp drop in the energy consistent with a first order transition 
is still present. The observed finite size transition temperatures agree well with those 




0.94 0.96 



Figure 6. The energy for Hduaiz on lattices ranging from 10 3 to 18 3 from left to right. 
The lines joining the data points are drawn to guide the eye. 

calculated by transforming the values from Fig. (4) using the decoration relation, e.g. 
for L = 10 we would expect (5 C = 0.5 x cosh -1 (exp(l. 275)) = 0.975, as found directly in 
the simulation. 

Further evidence for a first order transition with H duals, as noted above for the other 
dual Hamiltonians, can be garnered by looking at the energy histogram P(E) to discern a 
dual peak structure. In Fig. ^ close to the pseudocritical point for L = 10 at (3 = 0.972 
we can see clear evidence of two peaks. The Monte-Carlo simulations thus lend support 




-5.8 -5.6 -5.4 -5.2 -5 -4.8 




Figure 7. The energy histogram P(E) close to the pseudocritical point at /3 = 0.972 
on a 10 3 lattice. 



to the observation that H^uaiz and Hduan are related by a decoration transformation 
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through the agreement of the suitably transformed transition temperatures and also 
confirm the first order nature of the transition seen in H duals- There is no signal for 
the transition in the magnetic quantities (a), (r), (//) due to the gauge invariance and 
flip symmetries of the Hamiltonian, so in this respect too Hd ua i3 is similar to Hd ua i2 and 

Hduall- 

From the point of view of efficient simulation, an application of a principle of least 
effort suggests that the best adapted dual formulation for numerical work is probably 
that of the Ashkin- Teller like Hamiltonian of Hduaii in [5], since that has the minimum 
number of spin degrees of freedom to simulate. H dua i2 adds an additional flavour of 
vertex spin to this that is effectively a gauge degree of freedom and H dua i3 uncouples 
the four spin interactions in this with further edge spins. The physics of all three dual 
Hamiltonians is the same. 



6. Discussion 



To summarize, we have the following chain of equivalences between the various dual 
gonihedric Hamiltonians in 3d 

HdualS = ~Y1 ( a i U i0 a 3 + fkUijfij) ~ ( T i U ik T k + M^k) 

(ij) (ik) 

- Yl ( a 3 U !k°k + TjUf k T k ) 

(jk) 

— > (Un) Decoration — > 



H dua i2 = - Yj^^i^V-j ~ z2 TiTk ^ i ^ k ~ z2i (J i CTkT i Tk 

(ij) (ik) (jk) 

— > Gauge-Fixing — > (24) 



Hduaii = - - Y, T ' iTk ~ y^°"j°"fc r i r fc 

(ij) (ik) (jk) 

— > Non-Ising variables — > 



Hduaio = - 7", OiVj ~ Y TiTk ~ Y r] i T]k 

(ij) (ik) (jk) 

In the above we have listed the operations relating the various formulations in three 
dimensions. A variant of the classical decoration transformation in which edge spins are 
summed out relates H dua i3 to H dua i2- In transforming H dual3 -)■ H duai2 the coupling is 
therefore transformed as (3 = (1/2) In cosh(2/3). The gauge-invariant nature of H dua i3 
due to the presence of both edge and vertex spins leaves an echo in the vertex spin 
gauge symmetry of H dua i2, which in turn ensures the equivalence of H dua i2 and H dua ii 
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via a gauge-fixing. Allowing non-Ising spins gives a final equivalence between the dual 
models H dua n and H dua i and a standard duality transformation then takes us back to 
the original plaquette gonihedric Hamiltonian of H k= q. 

Gauge-invariant Hamiltonians such as H dua i3 have been employed in the past as 



models for open surfaces 13 , since the use of edge spins allows spin clusters to have 
free edges and seams [3). Indeed, the archetypal lattice gauge-matter theory, the 3d Z 2 
gauge-Higgs model fl4] 



H = -ft X>' r <^,) - & Yl UvUjkUkiUn , (25) 

(ij) [ijkl] 

has itself been used in such a context jl5J. In this paper the gauge spins U 1,2 ' 3 of H dua i 3 
are non-dynamical, so it would be an interesting extension of the current investigations 



to include a pure gauge term in the Hamiltonian in the manner of equ. (25) to observe 
the interplay between the anisotropic matter couplings and the gauge spins. 

Considering the original plaquette Hamiltonian H K=Q it is also possible to write 



down a gauge-Ising variant that allows open surfaces, which takes the form 16 

H = - /3 2 Y [{ a i U ij a j){ a kUklVl) + {ViUil(Tl){(JjU jk a k )} 
[ijkl] 

- f34j2(U ik U jk U M U u ) (26) 
[ijkl] 

where the matter couplings are dimer sums over the opposite edges of plaquettes. In 
this case the anisotropy seen in the dual Hamiltonians is not present. 

We close by repeating our observation in [6| regarding H dua i 2 and H dua n: it is a 
curious feature of these dual Hamiltonians, and as we have seen in this paper H dua i3 also, 
that they are all anisotropic in spite of being dual to an isotropic Hamiltonian. The flip 
symmetries that are a direct consequence of this, and the gauge symmetry manifest in 
H duals and inherited by H dua i 2 , play an important role in determining the properties of 
the models and how the different Hamiltonians are related. 

The coupling between the spin types and spatial directions in the dual gonihedric 
models is reminiscent of another class of anisotropically coupled Hamiltonians, the 



compass models, which exist in both classical and quantum forms 17 . In 3c? for instance, 



the quantum compass Hamiltonian is of the form 18 

H compass = ~J X ^ ^ °i ~ J V ~ J z ^i* 7 * ' ^ 

(ij) (ik) (jk) 

where the sums are again one dimensional and the a are now Pauli matrices. These too 
are know to display numerous interesting symmetries (191 , and have been found to have 



strong finite-size effects with periodic boundary conditions 20 . The parallels between 
these models and the gonihedric Hamiltonians merit further investigation, as do the 
potential numerical pitfalls involved in simulations of both. 
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